Can one measure nuclear matrix elements of neutrinoless double beta decay? 
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By making use of the isospin conservation by strong interaction, the Fermi Of /3/3 nuclear matrix 
element M F V is transformed to acquire the form of an energy-weighted double Fermi transition 
matrix element. This useful representation allows reconstruction of the total M% v provided a small 
isospin-breaking Fermi matrix element between the isobaric analog state in the intermediate nucleus 
and the ground state of the daughter nucleus could be measured, e.g., by charge- exchange reactions. 
Such a measurement could set a scale for the Of/3/? nuclear matrix elements and help to discriminate 
between the different nuclear structure models in which calculated M% v may differ by as much as a 
factor of 5 (that translates to about 20% difference in the total M ). 
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Neutrino is the only known spin-1/2 fermion which 
may be truly neutral, i.e., identical with its own antiparti- 
cle. In such a case one speaks about Majorana neutrino, 
to be contrasted with Dirac neutrino which is different 
from its antiparticle [HQ. Majorana neutrinos naturally 
appear in many extensions of the standard model (see, 
e.g., Q). Also, the smallness of neutrino masses (more 
than five orders of magnitude smaller than the electron 
mass) finds an elegant explanation within the see-saw 
model which assumes neutrinos to be Majorana parti- 
cles 01 . 

The fact that neutrinos have mass has firmly been es- 
tablished by neutrino oscillation experiments (for reviews 
see, e.g., Ref. 5]). However, the observed oscillations 
cannot in principle pin down the absolute scale of the 
neutrino masses. A study of nuclear neutrinoless double 
beta (Of/?/?) decay ^Elw — ► z ^,~El N - 2 + 2e~ offers a 
mean to probe the absolute neutrino masses at the level 
of tens of meV. 

Double beta decay is a rare decay process which may 
occur in the second order of weak interaction. It offers 
the only feasible way to test the charge-conjugation prop- 
erty of neutrinos. The existence of the Of/?/? decay would 
immediately prove neutrino to be a Majorana particle as 
assured by the Schechter-Valle theorem [||. The Of/?/? 
decay is strictly forbidden in the standard model of the 
electroweak interaction in which the lepton number is 
conserved, thus its observation would be of paramount 
importance for our understanding of particle physics be- 
yond the standard model @, d, U- 

The next generation of 0f/?/?-decay experiments 
(CUORE, GERDA, MAJORANA, SuperNEMO etc., 
see, e.g., Ref. [§] for a recent review) has a great discovery 
potential. Provided the corresponding decay lifetimes are 
accurately measured, knowledge of the relevant nuclear 
matrix elements (m.e.) M 0l/ will become indispensable 
to reliably deduce the effective Majorana mass from the 
experimental data. 

Two basic theoretical approaches are used to evalu- 



ate M 0v , the quasiparticle random phase approximation 
(QRPA) [13, U |, including its continuum version 
and the nuclear shell model (NSM) [l3(. There has been 
great progress in the calculations over the last five years, 
and now the QRPA Of/?/? nuclear m.e. of different groups 
seem to converge. However, the NSM M 0v are system- 
atically and substantially (up to a factor of 2) smaller 
than the corresponding QRPA ones. There is now an ac- 
tive discussion in literature on what could be the reason 
of such a discrepancy, a too small single-particle model 
space of the NSM or a neglect of complex nuclear config- 
urations within the QRPA. Even more striking is the dif- 
ference in the Fermi contribution to the total M 0l/ which 
can be up to a factor of 5 larger in the QRPA calculations 
than in the NSM ones. 

In view of this situation, it would be extremely im- 
portant to find a possibility to determine M 0l/ experi- 
mentally. There have been attempts to reconstruct the 
nuclear amplitude of two-neutrino /?/? decay (which ex- 
perimentally is very accurately known from the direct 
counting /?/?-decay experiments [9() from partial one- leg 
transition amplitudes to the intermediate 1 + states mea- 
sured in charge-exchange reactions [ljj]. However, such 
a procedure can consistently determine M 2l/ only if a 
transition via a single intermediate 1 + state dominates 
M 2v (the so-called single-state dominance). In the case 
of comparable contributions of several intermediate 1 + 
states the results from charge-exchange reactions cannot 
directly provide M 2v , since relative phases of the contri- 
butions cannot be measured. Pursuing the same way to 
reconstruct M 0u seems even more hopeless, since many 
intermediate states of all multipolarities (with a rather 
moderate contribution of the 1 + states) are virtually pop- 
ulated in the Of/?/? decay due to a large momentum of 
the virtual neutrino. 

The aim of this Rapid Communication is to suggest a 
way by which at least the Fermi component of M 0v can 
directly be measured, e.g., in charge- exchange reactions. 
For the derivation of the master expressions 61), © below 
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the well-known property of the Coulomb interaction to 
be the leadin g so urce of the isospin breaking in nuclei is 
exploited la, [l|| . Such a measurement of M F V could set 
a scale for the 0^/3/3 nuclear m.e. and help to discrimi- 
nate between different nuclear structure models in which 
calculated M F V may differ by as much as a factor of 5. 

We start our derivation by writing down the 0v(3f3 nu- 
clear m.e. in the closure approximation in which it ac- 
quires the form M 0u = (0 f \W 0u \0i) of the m.e. of a 
two-body scalar operator W au between the parent and 
daughter ground states |0j) and |0/), respectively. 1 The 
total O^/3/3-decay operator W 0v = g 2 A W^ T - g\rW F v is 
the sum of the Gamow- Teller and Fermi transition oper- 
ators 0: 



= ^2 Pv(r ab ) (g\(T a ■ (T b - 
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Here, the vector and axial vector coupling constants are 
gv — 1 and gA — 1.25, respectively, and P v (r ab = 
\r a — f%\) is the neutrino potential which in the simplest 
Coulomb approximation is just reciprocal of the distance 
between the nucleons: Pv(r ao ) — ^r— (for the sake of sim- 
plicity we have taken out the nuclear radius R from the 
usual definition of P v [7]). In this approximation 
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where T = r a is the isospin lowering operator, 



and V c = — 
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is the operator of 



Coulomb interaction between protons. Actually, only the 
isotensor component of the Coulomb interaction V c = 
e 2 T (2) 
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survives in the 



double commutator @. This isotensor Coulomb inter- 
action does contribute to the mean Coulomb field in the 
nucleus, but it is easy to see that any mean-field single- 
particle operator drops out of the double commutator ((2]) . 
Thus, the expression ([2]) is essentially determined by the 
residual (after separating out the mean-field contribu- 
tion) two-body isotensor Coulomb interaction. 

The total nonrclativistic nuclear Hamiltonian H to t 
consists of the total kinetic energy of nucleons and 
the strong and Coulomb two-body interactions between 



1 Using closure of the states of the intermediate nucleus z+ ^Eljv— l 
which are virtually excited in /3/3-decay would be an exact proce- 
dure if there were no energy dependence in the 0v[3(3 transition 
operator. A weak energy dependence of the operator leads in 
reality to a "beyond-closure" correction to the total M " of less 
than 10%. 



them: H tot = T + H str 
actly isospin-symmetric 
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(we shall quan- 



tify later the accuracy of this assertion but it is well 
known that the isospin-breaking terms in H str are in fact 
fairly small [H,[ll]), one has 
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Here, the sum runs over all + states of the intermedi- 
ate (N — 1,Z+ 1) isobaric nucleus, Q s = E s — (E 0i + 
Eo f )/2 represents the excitation energy of the s'th in- 
termediate + state relative to the mean energy of the 
ground states of the initial and final nuclei. To ac- 
count for the isospin-breaking part of H s t r , SMp 11 = 



3 <0/l 
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\0i) should be subtracted from 

r.h.s. of Eq. {4}. 

Among all the intermediate + states, the isobaric 
analog state (IAS) dominates the sum J3J). In fact, 
(IAS\f-\0i) w - Z is the largest first-leg transition 
m.e. (a few percents of the total Fermi strength N — Z 
may go to the highly-excited isovector monopole reso- 
nance (IVMR) since the IAS and IVMR get mixed mainly 
by the Coulomb mean field). Similarly, the second- leg 
Fermi transition dominantly populates the double IAS 
(DIAS) in the final nucleus. Due to the isotensor part of 
the Coulomb interaction (which also gives the only con- 
tribution to the double commutator ([2|)), the final g.s. 
gets an admixture of the DIAS where the correspond- 
ing mixing m.e. is (OADIAS) = - ^f\ V c\ DIAS ) ^ witn 

Edias 

E DIAS w 2cu IAS . Thereby, one gets (0 f \f-\IAS) ± 0. 

Other quantitative arguments for the dominance of the 
IAS in the sum ([4} follow from the representation of the 
double commutator: 



T- [f- V£] 



{f-)\(f-j 



% - 2T~V£f- 



It is clear that the first term V c 
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(T ) dominates the 
since the other m.e., be- 



(0/1 

cause of T + |0/) « (with a small deviation from zero 
originating from an isospin symmetry violation effect, 
caused mainly by the Coulomb mean field), contain ad- 
ditional suppression as compared with the leading term 

(0/lVg (f-f \0i) = (0 f \V c \DIAS)(DIAS\ (f-f |0*>. 

Thus, M F V is determined by the amplitude of the dou- 
ble Fermi transition via the IAS in the intermediate nu- 
cleus into the ground state of the final nucleus where 
the second Fermi transition amplitude is due to an ad- 
mixture of the DIAS in the final nucleus to the ground 
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state of the parent nucleus: (Q f \f-\IAS){IAS\f-\Oi) = 
(O f \DIAS)(DIAS\f-\IAS)(IAS\f-\0 % ). Finally, one 



can write 



■u>iAs(O f \T-\IAS)(IAS\f-\Oi 
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Therefore, the total M F V can be reconstructed accord- 
ing to Eq. |(SJ|, if one is able to measure the AT = 2 
isospin-forbidden m.e. (Of\T~\IAS), for instance in 
charge-exchange reactions of the (n,p)-type (also the 
same m.e. determines Mp v , but it would be much more 
difHcult to extract it}. Using the QRPA calculation re- 
sults for M F V [l(J, [IH, this m.e. can roughly be esti- 
mated as (Of\T~\IAS) ~ 0.005, i.e. about a thousand 
times smaller than the first-leg m.e. (IAS\T~\0i) s» 
y/N - Z. This strong suppression of (0 f \f~\IAS) re- 
flects the smallness of the isospin violation in nuclei. The 
IAS has been observed as a prominent and extremely nar- 
row resonance and its various features have well been 
studied by means of (p,n), ( 3 He,t) and other charge- 
exchange reactions, see, e.g. This gives us hope 
that a measurement of (0f\T~\IAS) in the (n,p) charge- 
exchange channel might be possible. More generally, 
a measurement by whichever experimental mean of the 
AT = 2 admixture of the DIAS in the final ground state 
would be enough to determine M F . 

A qualitative analysis of the physics involved in cal- 
culations of M F V can be conducted further. One can 
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define an operator 
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taincd by the substitution of — by a constant ■= 

Tab _ R 

in the definition of the isotensor Coulomb interaction. 
Such an operator Vq is diagonal in the basis of isospin 
eigenstates and does not mix in the first order the 
DIAS and the final ground state, (DIAS\%\0 

0. The matrix element — (0/| f~ 
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a large factor 



— — <C 1 smaller than the absolute value of the r.h.s 
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in Eq. (HI). Thus, by subtracting from V$ in Eq. © 
only a small change in M F V is introduced. Such a subtrac- 
tion with an appropriate choice of R ~ R allows to cut 
off the contribution to M F U from the long internucleon 
distances, where — has a smooth behavior and which 
are relevant for the Coulomb mean field. Therefore, the 
major contribution to M F V should come from the short 
distances where the gradient of — is the largest. This 
provides a natural qualitative explanation of the numer- 
ical results of both the QRPA and NSM [HI El which 
consistently show the short-range character of the partial 
r-dependent contribution to M 0ly . 

Of course, by measuring only M F V one does not get 



the total m.e. M 0l/ but rather its subleading contribu- 
tion. However, knowledge of M F V itself brings a very 
important piece of information. For instance, it will 
allow to investigate the A dependence of M F V . Also, 
it can help a lot to discriminate between different nu- 
clear structure models in which calculated M F V may dif- 
fer by as much as a factor of 5. In addition, the ratio 
Mp? /Mfvp may be more reliably calculable in different 
models than M F V and Mgy separately. Let us put for- 
ward here some simple arguments in support of the lat- 
ter statement. Since only small internucleon distances 
determine M 0v , then only nucleon pairs in the spatial 
relative s-wave must dominantly contribute to the m.e.. 
The isotensor Coulomb interaction only couples T = 1 
pairs which must then be in the state with the total 
spin S = to assure antisymmetry of the total two- 
body wave function. Because of this and the fact that 
(7i ■ eralS 1 = 0,T = 1) = -3|5 = 0,T = 1), a natural 
estimate for the Gamow- Teller m.e. is M^ T = -ZM F V 
provided the neutrino potential is the same in both F 
and GT cases. The high-order terms of the nucleon weak 
current which are present in the case of the GT m.e., but 
absent in the F m.e., change a bit this simple estimate 
to M$y T /M F v w -2.5. Also, an uncertainty of few per- 
cent may come from the difference in the mean nuclear 
excitation energies in the F and GT cases. It is worth 

w 
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noting that the recent QRPA results 
good correspondence with these simple estimates. 

Here, we want to estimate possible corrections to the 
simplest closure approximation discussed above. Due to 
universality and conservation of the vector current, all 
the corrections of the vector current vertices should be 
the same independently of which virtual particle, neu- 
trino or photon, is exchanged between them. This is true 
for the effects of short-range correlations and the finite 
nucleon size. A small difference of a few procent in the 
realistic potentials may arise from different mean nuclear 
excitation energies while exchanging the neutrino or pho- 
ton but this effect seems to be rather reliably calculable. 
Another difference can arise from those corrections to the 
propagator of the virtual photon, as for instance the vac- 
uum polarization correction, that are missing in the case 
of the virtual neutrino. The effect of the the vacuum po- 
larization is about 0.5% and can simply be accounted for 
by a proper renormalization of the electron charge. 

The effect of isospin nonconservation in the strong two- 
body interaction can be estimated to be at the level of 

2%-3% mm. One can then directly compare the ra- 
dial dependencies of the isospin-breaking part of the two- 
body strong interaction in the S = 0, T — 1 channel and 
the Coulomb interaction within the relevant short range 
of 1-2 fm to find the dominating source of the isospin 
breaking. Following Ref. [19j], one can approximate the 
radial dependence of the isospin-breaking strong two- 



body central potential as (0.02-0.03)x^ ; 



{h = c = 
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1). With 
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0.08 one arrives at the conclusion that this 



source of the isospin non-conservation must be about 20- 
30 % of that caused by the Coulomb interaction. Though 
there are rather large relative uncertainties in calculating 
the isospin-breaking part of the two-body strong interac- 
tion, by assuming that this correction could in principle 
be evaluated with a moderate accuracy of 30 %, a resid- 
ual uncertainty of only 10 % in Mp" is thereby induced. 

Thus, the main message of this Rapid Communication 
that, at least in principle, M F V is measurable remains 
intact in the most realistic situation (though minor cor- 
rections may be needed). 

To conclude, we have shown in this Rapid Communi- 
cation that the Fermi 0i//3(3 nuclear m.e. can be recon- 
structed if one is able to measure the isospin-forbidden 
Fermi m.e. between the ground state of the final nu- 
cleus and the isobaric analog state in the intermediate 
nucleus, for instance by means of charge- exchange reac- 
tions of the (ra,p)-type. Knowledge of M F U would bring 
a quite important piece of information on the total Qv(3f3 
nuclear m.e.. Simple arguments show that the estimate 
M^V T /M F V sa -2.5 should hold. Also, such a measure- 
ment can help to discriminate between different nuclear 
structure models in which calculated M F V may differ by 
as much as a factor of 5. 
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